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ON THE DISTRIBUTION OF THE GREATEST COMMON 
DIVISOR OF GAUSSIAN INTEGERS 

TAI-DANAE BRADLEY, YIN CHOI CHENG, AND YAN FEI LUO 


Abstract. For a pair of random Gaussian integers chosen uniformly and in¬ 
dependently from the set of Gaussian integers of norm x or less as x goes to 
infinity, we find asymptotics for the average norm of their greatest common 
divisor, with explicit error terms. We also present results for higher moments 
along with computational data which support the results for the second, third, 
fourth, and fifth moments. The analogous question for integers is studied by 
Diaconis and Erdos. 


1. Introduction 

In this paper, we study questions related to the size of the greatest common 
divisor of pairs of randomly chosen Gaussian integers. In particular, in Theorem [TJ 
we first calculate the probability that a pair of random Gaussian integers, chosen 
uniformly and independently from the set of all Gaussian integers with norm x or 
less, has greatest common divisor ±k or Hk for a fixed Gaussian integer k. The 
main term for this probability in the case where k =1 was first given by Collins and 
Johnson |3} Theorem 8]. We refine their results by providing the expression for 
the more general case in addition to giving an explicit error term for all cases. In 
Theorem [2] we derive the expected norm of the greatest common divisor between a 
pair of Gaussian integers with norm x or less. Finally, in Theorem [3] and Conjecture 
[TJ we present an expression for higher moments of the norm of the greatest common 
divisor between a pair of Gaussian integers with norm x or less. We expect our 
results to generalize to principal ideal domains without too much difficulty. More 
generally, our results should hold for the ring of integers in an algebraic number field, 
though our techniques will need to be modified to deal with class number greater 
than one and infinite unit group. We expect the ideas in the recent preprint of [8] 
could help address this question and would be an interesting direction to explore 
further. Of further interest are function field analogues. Some interesting results 
in this direction may be found in [9]. 

Similar questions have also been studied for the case of rational integers. Orig¬ 
inally, Mertens [7 proved in 1874 that the probability that a pair of rational inte¬ 
gers chosen uniformly and independently at random from { 1 , 2 ,..., a;} are relatively 
prime is asymptotic to l/£(2), as x tends to infinity, where £ is the Riemann zeta 
function. In 1956, Christopher Theorem 1] generalized Mertens’ result by finding 
the probability that two integers have greatest common divisor k for a fixed k larger 
than 1. An asymptotic expression for the moments of the greatest common divisor 
was first derived by Cesaro in 1885 [T], and Diaconis and Erdos later extended his 
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work by explicitly calculating the error term [ 4 ], Theorem 2 ]. In particular, the 
expected value for the greatest common divisor between a pair of random integers 
chosen independently and uniformly from the set {1,2,..., a:} is 

(!) logx + 0(1) 

while the n th moment is given by 

(2) Ln{<fri)“ 1 } +0(l ”^ logl) for " a2 ' 

The goal of the present paper is to show that equation ([I]) has an analogous 
counterpart in the ring of Gaussian integers as stated in Theorem [ 2 ] at the end 
of this section. Further, we show that equation | 2 | also has an analogous form 
as presented in Theorem [ 3 ] and Conjecture [T] Before proceeding, we first give the 
following preliminary definitions and remark. 

Definition 1. The norm of a Gaussian integer a = a + bi for rational integers a 
and b is defined by N(a) = a 2 + b 2 . 

Most of our results will be in terms of the norms of Gaussian integers and not 
the integers themselves. 


Remark 1 . Given two Gaussian integers 77 and 71, a greatest common divisor, de¬ 
noted (77, 77), is defined to be a Gaussian integer k such that k is a divisor of both 77 
and 77, and if there is any other common factor between 77 and 77, then it must also 
be a factor of k. From this definition, it becomes clear that (77, 77) is unique only up 
to its associates. In other words, (77, /i) = n, —k, in, and — in. Our calculations, 
however, will be performed via ideals for reasons that will soon become apparent. 
For a Gaussian integer 77, ti is the ideal such that n = (77) = (—77) = (777) = (—777), 
and the norm of n is defined by IV(n) = N{rj). Accordingly, the definition of the 
greatest common divisor for a pair of ideals is: 


Definition 2 (Greatest Common Divisor of Two Ideals). For a ring R , let n, m C R 
be ideals. The greatest common divisor (n, m) is defined to be the ideal .ft C R which 
satisfies the following: 

(1) n C A and m C ft 

( 2 ) if there exists some ideal a C R such that n C 0 and m C a, then ft C a. 

In other words, (n, m) is the smallest ideal that contains all the elements of both 
n and m. When applied to the ring of Gaussian integers, a Dedekind domain, it is 
clear that (n, m) is unique. 


Definition 3 (The Dedekind Zeta Function of Q(i)). For the number field Q(7), 
the complex-valued Dedekind zeta function is defined for Re(s) > 1 by 


«( s ) = H 


aC Z[i] 


iV(a) 6 


(a 2 

(a,&)ez v 

(o,6)#(0,0) 


b 2 ) s 


where the first summation is over the nonzero ideals a of the ring of Gaussian 
integers Z[i]. 

In order to find the expression for the expected norm of a greatest common 
divisor between a pair of Gaussian integers of norm x or less, we will first derive 
the necessary probability distribution function of Theorem [T] 
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Theorem 1. Let n and m be nonzero ideals chosen independently and uniformly 
at random from the set of ideals in Z[z] with norm x or less. The probability that 
(n, m) = M. is 

+ ° (a; 2 /37V(j?) 4 / 3 ) ' 

This probability will allow us to calculate the expected norm of the greatest 
common divisor between a pair of ideals: 


Theorem 2. Let n and m be nonzero ideals chosen independently and uniformly 
at random from the set of ideals in Z[z] with norm x or less. The expected norm of 
the greatest common divisor of n and m is 

k irn' mx+0{1) - 

We will then prove the following result regarding the n th moment for n > 2 : 


Theorem 3. Let n and m be nonzero ideals chosen independently and uniformly at 
random from the set of ideals in Z[z] with norm x or less. For n > 2, there exists 
a constant c n £ ffi. such that 

E x {N(n, m)"} ~ c n ai n_1 

where E x {N(n, m) n } denotes the n th moment of the norm of the greatest common 
divisor of n and m. 


Lastly, we will present numerical data which provide strong evidence for the 
following conjecture regarding the constant of Theorem [ 3 ] for all n > 2 : 


Conjecture 1. For n> 2, 

E x {N(n, m)"} 


2Cq(») ( n ) 
n(n + 1) 1 Cq (i)(n + 1) 


- l>x 


n — 1 


The proof of Theorem [l] will be given in Section [ 2 ] and that of Theorem [ 2 ] will 
be given in Section [ 3 ] Finally, in Section [ 4 j we prove Theorem [ 3 ] and present 
Conjecture [I] along with computational data which support the conjecture for the 
second, third, fourth, and fifth moments. 


2. Probability Distribution Function 

Before deriving the expression for the probability of Theorem [TJ we first define 
the following two functions: 

Definition 4 (The Mobius Function). For an ideal n, the Mdbius Function ^i(n) is 
defined by: 

fl, if n = (1); 

/ i ( n ) = \ (—!)*, if R = P1P2 • • • p t f° r distinct prime ideals p, : ; 

[0, otherwise. 

We will use the following identity 

= 

D|n 


( 3 ) 


1 , if n = (1); 
0 , if n 7^(1) 
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as well as the generating function 

Mn) = 1 

N (n) s C Q>(i){s) 


E 

n CZ[i\ 


for Re(s) >1. 


Definition 5 (The Sum of Two Squares Function). For n € Z, let the sum of two 
squares function r(n, 2) represent the number of ways that n can be expressed as a 
sum of two squares. Thus, 

r(n, 2) = i#{a C Z[i\ : N( a) = n}. 

We will need the result of Sierpinski (11 (for a statement in English see flOj equation 

( 1 )]) 

( 4 ) 


E r ( n > 2 ) = 


7 TX - 


0(xz 


n—l 


The error term 0(a; 1 / 3 ) has been improved by Huxley 6j to 0(x™6 +e ) } but the 
former is sufficient for our purposes. We shall also use the following m 


( 5 ) 


Y r ( n? 2) = 7 r(5 + logx) + 0{x 1//2 ) 

n—l 


where S denotes Sierpinski’s constant S ~ 
expressions 0 P- 123] 

S = — lim (4 ((z)(3(z) — 


2.58/7T. This also has the alternate 

^(l) 


z-1 


= 7 + 


m 


where f3(z) is the Dirichlet beta function and 7 is the Euler-Mascheroni constant. 

With these functions at hand, we may now proceed to calculate the desired 
probability. To do so, we will need two preliminary results. The first is the total 
number of pairs of ideals generated by Gaussian integers with norm at most x. The 
second result is the number of those pairs which have greatest common divisor H. 
The expressions for each of these are derived in the following two lemmas. 

Lemma 4. The total number of pairs of nonzero ideals n and m in Z[i\ with norm 
x or less is 


:k 2 x 2 

16 


+ 0(x A/3 ). 


Proof. Let n and m be nonzero ideals. Then 

#{n,m C Z[i\ 2 : IV(n), JV(m) < x} = ^ E 1; 


nCZ[i] mCZ[i] 

TV (n ) < x lV(m)<x 


and we may rewrite this as 


1 

16 


L*J 


[x\ 


E r ( iV ( n )> 2 ) E r ( 7V ( m )^ 2 ) 


IV(n)=l 

which by equation @ equals 


AT(m)=l 


1 

16 L 


TTX + 0(lzj 1 ^ 3 ) 
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Further, since 0 (|_ x\) = 0(x) we may expand \nx + 0( a; 1 / 3 )] 2 and obtain ir 2 x 2 + 
2irx ■ 0(x 1 / 3 ) + 0(x 2 / 3 ) which reduces to n 2 x 2 + 0(x 4 / 3 ). Thus, the total number 
of n and m with norm at most x is 

~-2 ™2 

+ o(- 4/3 ). 

□ 


Lemma 5. The total number of pairs of nonzero ideals n and m in Z[i\ with norm 
x or less having greatest common divisor A is 

ir 2 x 2 ( a; 4 / 3 \ 

i 6 Co(i )( 2 ) fc2 + U 4/ v 

where k = N(&). 

Proof. Let n and m be nonzero ideals. The number of pairs of n and m with norm 
x or less which are relatively prime is 

#{n, m C Z[i] 2 : N(n),N(m) < x and (n, m) = ( 1 )} 

= E Ei 

nCZ[i] mCZ[i] 

N (n) <.x iV(m)<a; 

(n,m) = (l) 


= E E E vi*) 

nC Z[i] mCZ[i] t)CZ[i] 

N(n)<x N (m)<x 0|(n,m) 

where in the last line we have used identity ©>■ Reindexing with n = On' and 
m = Om' where the norms of n' and m' range from 1 to x/N(d), we may rewrite 
this as 

E E E 1 

t)CZ[i] n'C m'C Z[i\ 


-t L J >> ) J J 

= 16 E E r ( iV ( n, ),2) E r ( jV ( m ') I 2). 

DCZ[i] Nn' = l Nxn'=l 

AT(D)<x 


As in Lemma [21 this reduces to 


1 

16 


E 

t>C Z[i] 

N( 0)<x 


N(vy 


o 



We then distribute the summation to obtain 















6 


TAI-DANAE BRADLEY, YIN CHOI CHENG, AND YAN FEI LUO 


To evaluate the main term, we call on the generating function Dn C z[»] N(n )° = 
1/Cq(j)( s ) for Re(s) > 1, to see that 

MM _ i _ y' y" MM 

Z—* N Ci))" 2 CorMM Z—* -NYd) 2 

ocz[i] v syojv ) n=x +i ac z[i] v ' 

JV(0)<3: N(Q)=n 


which implies 
1 


Cq(»)( 2 ) 


- E 

DCZ[i] 

N(i3)<x 


MM 


N(py 


OO -J 

i E?E> 

n—x +1 DCZ[i] 

N(a)=n 


1 

4 


-(n, 2) 


n—tc+1 


Now we note that r(n, 2) < 4cr 0 (M = o(n e ) for all e > 0, where a 0 represents the 
number of divisors of n. Thus \ Y^™= x +i - E^U+i = oix^ 1 ) and so 


1 


- E 


MM 


< o(x e ~ 1 ) 


4 Z-^n=x-\-1 n 2 
e—l\ 


^)( 2 ) 0 ^(M 2 

iV(0)<x 

For the error term of ([6]), we have 
/ X \ 4/3 


or 


E 

0CZ[j] 

N(i>)<x 


MM = i 

^(M 2 Cq W (2) 




E 

DCZ[i] 

iV(D)<x 


U(M 


Y' 1 y" i — - V r ^ n ’ ^) 

— Z-j n 4/3 2^ 4 L n 4/3 

n=l OCZfzl n=l 


«(o)= 


and again use the bound r(n, 2) < o(n e ) to see that j L)n =1 ^ 4 / 3 ^ < En=i 
which equals o(x e_1 / 3 ) + o(l). From this it is clear that O ^x 4//3 J2n=i 
0(o(x 4 / 3 )) = 0(x 4 / 3 ). Thus ^ becomes 


o(n e 4 / 3 ) 

r(n,2) \ _ 

n 4 / 3 y 


16Cqw(2) 


+ ofs*" 1 ) + 0(x 4 / 3 ) 


which allows us to conclude 

2 2 

#{n, m c Z[i] 2 : iV(n), AT(m) < x and (n, m) = (1)} = J X , . + 0(x 4/3 ). 

lb CQ(M 2 ) 

Having counted the number of relatively prime n and m within a given norm, we 
can now reindex to obtain the number of them which have (n, m) = J?. Letting 
n = n'j? and m = m'J?, we see that n' and m' are relatively prime if and only if n 
and m have J? as their greatest common divisor. Hence, the number of relatively 
prime pairs n' and m' with norm y or less must be equivalent to the number of pairs 
n and m, with norm yk or less (where k = having greatest common divisor 

&. Thus, 


#{n, m C Z[«] 2 : iV(n), N(m) < x and (n, m) = £} 

= #{n', m' C Z [i] 2 : -/V(n'), iV(m') < j and (n', m') = (1)} 

rC 

n 2 x 2 /x 4 / 3 \ 

~ 16C Q(i )(2)fc2 +C; l vfc 4/3j- 

□ 
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At last, the probability that n and m, having norm at most x, will have greatest 
common divisor A is defined to be the number of pairs of ideals of norm x or less 
which have greatest common divisor A divided by the total number of pairs of ideals 
of norm x or less. Thus, by Lemmas [4] and [5j 

(7) 7A{n, m c Z[i] 2 : N(n),N(m) < x and (n, m) = U} 

A2 _ +o (—V 

16 ' { J 16Coru(2)fc 2+ U 4/ vJ' 


We can rewrite 


-l 


_^t + 0(x 4 / 3 ) 

to 167r _2 a: _2 [l + 0(a; -2 / 3 )] since [1 + f(x)] 
towards 0 as x approaches infinity. 

Line ([TJ) then becomes 


7T 2 x 2 1 + 0{x 2/3 ) 


1 + 0(x~ 2 / 3 ) 

1 


16 Cq «)(2 )k 2 

as 167t _2 .t _2 [l + 0(a; -2 / 3 )] 1 which is equal 
= 1 + 0(f(x )) for /( x) tending 

fx 4 / :i \ 


n 2 x 2 


CQ(i)(2)fc 2 


O 


U 4/3 ) 


1 


Cq(»)( 2 )fc ; 


+ o 


Cdj(i)(2)fc 2 

1 


o 


1 


a; 2 / 3 /c 4 / 3 


O 


a ; 2 / 3 fc 4 / 3 
\x 2 / 3 k 2 


O 


1 


x 4 / 3 fc 4 / 3 


or finally 


+ O 


( 1 


C<Q(i)(2)A 2 ' ~ \x 2 / 3 k 4 / 3 / 

completing the proof of Theorem [I] The following corollary is a direct consequence 
of Theorem [l] for the special case when A = (1). 


Corollary 6. The probability that a pair of Gaussian integers with norm x or less 
are relatively prime is 

Cq^2) +0 (^)‘ 

In effect, Corollary [6] tells us that for x large, the probability that two Gaussian 
integers are relatively prime is asymptotic to [CQ(i)(2)] _1 as x tends towards infinity. 
This is in agreement with the work of Collins and Johnson who state the probability 
as [Cq W (2)]- 1 = [£(2)L(2, x)] _1 « 0.6637, where L(2,%) is a Dirichlet L-series and 
X the primitive Dirichlet character modulo 4. 


3. Expected Value 


Having derived the probability distribution function found in Theorem [T] we 
are ready to find an expression for the expected value of our random variable, 
N( n, m) = k, where the norm of n and m ranges from 1 to x. To do this, we must 
express our probability in terms of k as well. The modification is simple, however. 
Since the number of ideals with norm k is equivalent to r(fc, 2)/4, the probability 
that the greatest common divisor of n and m has norm k must be 


P x {iV(n,m) = k} 


2) / r(k, 2) \ 

4c qw (2)p + U 2/3 fc 4/3 y' 
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Then, by definition of expected value 

X 

E X {N( n, m)} = ^ k ■ P X {N{ n, m) = k} 

k=1 


= E fc 


k=1 


( 8 ) 


r (k, 2) + Q ( r(k,2 ) ^ 

Q(i '"" n 

X 

E 


.4Co(i)(2)fc 2 

r(fc,2) 


\ a; 2 / 3 fc 4 / 3 / 

1 


+ O 


E 


r(/c,2) 


2-2/3 fcl/3 

fe=l 


4 C<o(i)( 2 ) fe=1 

Using Stieltjes integration by parts to evaluate the error term, we obtain 


E ! ^7? = :r 1/3 E r ( fc ’ 2 )- 4 -/ ( 7rfc+0 ( fcl/3 )) (-^ fc 4/3 ) * 


fc=i 


= —x 2 / 3 + O(logx) 

which implies 0(x -2 / 3 ]Cfc=i ) = 0(1 + x -2 / 3 logx) = 0(1). The main term 
of ([8]) can be rewritten using Sierpinski’s identity from equation <©• Thus the 
expected value is equal to 


or 


1 

4 Cnj(i) (2) 


7 t(S + logx) + O 



7 r 

4 Cu(i)(2) 


logx + 0(1). 


This completes the proof of Theorem [2] 


+ 0 ( 1 ) 


4. Higher Moments 


At last, we show that there exists some constant c n £ R. such that the main 
term of the n th moment of iV(n, m) must be of the form c n x n_1 for n > 2. Let 


N(n),N(m) < x with (n, m) 


= 8. and restrict N(&) to the interval 


X X 

3 + 1 ’ 3 


. We 


may then write n = n'A and m = m'A where (n',m') = (1). 
on the norm of M allows us to see that N(n'),N(m') < x ■ 
N(W), N(m') < j. Now define 


The restriction 
■ which implies 


fU) = #{(n',m') C 1[i} 2 : N(n'),N(m') < j and (n',m') = (1)} 


for j £ N. By Lemma [5] 


f(j) 


1fi r 3 (0] + 0(i 4/3 ) 

16Cq(»)( 2 ) 

0(J 2 )- 


Our reindexing above shows that this expression for f(j) also gives us the number 
of pairs of ideals with norm x or less having greatest common divisor where 
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< N(ft) < j. Thus the n th moment of N (n, m) is given by 


(9) E X {N( n,m)"} = 


[7r 2 a; 2 /16 + O(o; 4 / 3 )] 


E /O') E n &y 

3 =1 &CZ[i] 

*<N(Si)<Z 
j + 1 v J 


We next turn our attention to the inner sum of ([9]). First note that 


E N(ft) n = - E k n r(k, 2) 


LfJ 


fiCZ[i] 


fc =rjfii 

where k = N(ft). Then Stieltjes integration by parts yields 

LfJ LfJ \M 


fe=rjfrJ 


n+1 


4(n +1)‘ 


jn +1 Q- _|_ ]^ra+l 


o 


E k n r(k,2) = \x/j\ nS ^r(k,2) - E r (M) ~ ^ 1 (?rt + 0(t 1/3 ) 

af 1 fc—1 -l 4 ” fc=l -'rfprJ 

+ 1/3 


The numerator of E x {N(n, m) n } is now equal to 


( 10 ) 


-x 


n+l 


E°o 2 ) 


4(n+l) 

The sum on the left is 


(j + l) n+1 ~ j n+1 

jn+l^j _j_ \)n+l 


X 


n+1/3 


E°(j 2 )° 


E°o- 2 )°f 


i 


j=i 


Vi(i + i ) n+1 


/=i 


i=i 

1 


n+l/3 


^°Vi- 1 (i + i) n+1 


which is bounded above by ff) 0(1/j n ). For a; tending toward infinity and n > 2, 
this converges to some constant c4 £ R. A similar argument shows that the second 
sum of ( fl0| is likewise convergent for n > 2. We thus conclude that the main term 
of .E x {iV(n, m)™} is of the form c' n x n+1 . 

At last, we divide this by the total number of pairs of ideals n, m with norm at 
most x to obtain the main term of the n th moment of iV(n, m) for n > 2: 

J rr n +1 i 

_ L n x _ _ „ x n- 1___ 

[7r 2 a; 2 /i6 + 0(a; 4 / 3 )] " 1 + 0(x~ 2 / 3 ) 

where c n = lQc' n /n 2 . Since [1 + 0(x~ 2 / 3 )]~ 1 = [1 + 0(x ~ 2 / 3 )], it follows that for x 
tending to infinity 

E x {N(n,m) n } ~ CnX n -\ 

With this, we bring the proof of Theorem [3] to an end and close by restating our 
conjecture regarding the constant of E x {N(n, m)”} for all n > 2. We also include 
numerical evidence below which provides support for the conjecture in the cases 
when n = 2,3,4 and 5. 
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Conjecture 1. For n> 2, 

E X {N{ n,m) n } 


2 Cq (i){n) 

7 T (n + 1) 1 + 1) 


— 1 i 


„n— 1 


Using Matlab, we first compiled a list of all pairs of Gaussian integers in the first 
quadrant with norm a or less and used the Euclidean Algorithm to find all possible 
greatest common divisors. We determined the n th moment by raising the norm of 
each greatest common divisor to the n th power, summed the terms together, and 
then divided the result by the total number of pairs of Gaussian integers in the first 
quadrant with norm x or less. We have graphed the results in Figures [I]-[4] below 
for the cases when n = 2, 3,4 and 5 with x = 50, 000. In Table[l] we have listed the 
main term of the best fit curve corresponding to each graph as compared against 
the conjectured main term for each value of n. 


Table 1. The main term of the n th moment of the norm of the 
greatest common divisor of pairs of Gaussian integers with norm 
at most x. 


Moment (n) 

Numerically Derived Term 

Conjectured Term 

2 

0.63952a; 

0.67364a; 

3 

0.37018a; 2 

0.37444a: 2 

4 

0.27238a: 3 

0.27309a 3 

5 

0.21914a; 4 

0.21928a 4 


Figure 1. The graph of E x {N(n,m) 2 } for 1 < x < 50,000. The 
best fit curve is 0.63952a+ 0.5753. 
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Figure 2. The graph of E X {N( n, m) 3 } for 1 < x < 50,000. The 
best fit curve is 0.37018a; 2 + 0.69337a; - 584.8498. 



Figure 3. The graph of E x {N(n, m) 4 } for 1 < x < 50,000. The 
best fit curve is 0.27238a; 3 +0.80149a; 2 

3723.1433a;+12324561.4508. 
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Figure 4. The graph of E X {N( n, m) 5 } for 1 < x < 50,000. 
The best fit curve is 0.21914a; 4 + 0.92436a; 3 — 9773.8223a; 2 + 
92150266.2382a; - 190551355734.3794. 

,J7 
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